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Abstract 

Exchange interaction between conduction electrons and magnetic moments at magnetic interfaces 
leads to mutual conversion between spin current and magnon current. We introduce a concept of 
spin convertance which quantitatively measures magnon current induced by spin accumulation 
and spin current created by magnon accumulation at a magnetic interface. We predict several 
phenomena on charge and spin drag across a magnetic insulator spacer for a few layered structures. 

PACS numbers: 72.25.Mk, 75.30.Ds 



I. INTRODUCTION 



In spintronics, spin current, which is conventionally defined as the difference of electric 
currents of spin-up and spin-down conduction electrons, plays a pivotal role in propagating 
spin information from one place to another. Many spin dependent properties, such as giant 
magnetoresistance [U [2] , spin transfer torques [3j H] and spin Hall effect j5j [6] , are directly 
related to spin current. Spin current has several unique properties compared to charge 
current: 1) it is considered as a flow of angular momentum while the conventional current 
is a flow of charge, 2) the total spin current is not a conserved quantity even in the steady 
state condition; it can be transferred and/or lost due to spin-dependent scattering, and 3) 
spin current has both transverse and longitudinal components whose decaying length scales 
are quite different in a ferromagnetic medium. Recently, the concept of spin current has 
been extended to spin wave current since spin waves carry angular momenta as well [7j. 
There are two types of spin wave currents. One is magnetostatic wave propagation [7HTT] 
for which the classical magnetization is temporal and spatially dependent. An example is 
a moving domain wall driven by a magnetic field or by an electric current. Although such 
magnetostatic spin waves may carry angular momentum, they are not quasi-particles in that 
there are no particle numbers associated with these waves. The other spin wave current is 
a true quasi-particle current known as magnon current. A magnon is a quantum object 
(particle) that represents low excitation state of ferromagnets. In equilibrium, the number 
of magnons AT q can be cast into a simple Boson distribution N® = [e E ^ kBT — where 
is the magnon energy. Similar to the electron spin, each magnon carries an angular 
momentum —Ti. In thermal equilibrium, there is no magnon current since there are an equal 
number of magnons moving in all directions. 

In our earlier paper [T2] . we showed that the non-equilibrium magnon accumulation and 
magnon current can be treated semiclassically, similar to the spin transport properties of 
conduction electrons. We found that the non-equilibrium electron spin current in metal 
can convert into a magnon current of a magnetic insulator through the interfacial exchange 
interaction. The magnon current then subsequently diffuses inside the magnetic insulator. 
The magnon diffusion process may be described by the diffusion equation. Among other 
things, we predicted that an electric current applied in one metallic layer can induce an 
electric current in another metallic layer separated by a magnetic insulator via magnon 
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mediated angular momentum transfer. In this paper, we extend our theory to include a 
general boundary condition for the spin convertibility at metal | magnetic-insulator interfaces 
and then, we calculate the electric drag in a few realizations. The paper is organized as 
follows. In Sec. II, we summarize the general boundary conditions at the interfaces between 
metals and magnetic insulators. In particular, we introduce a quantity, named as spin 
convertance, which quantitatively characterizes conversion effectiveness among spin/magnon 
accumulation and magnon/spin current at a magnetic interface. In Sec. Ill, we calculate the 
spin convertance by using the microscopic s-d exchange interaction. In Sec. IV, we present 
the general solutions for several layered structures with a magnetic insulator layer (MIL) 
and discuss some limiting cases. Finally, we summarize our results. 



II. SUMMARY OF BOUNDARY CONDITIONS 



We consider a simple bilayer consisting of a metallic layer in contact with a MIL. The 
angular momentum in a metal is carried by conduction electrons while in a MIL, it is carried 
by magnons. In the semiclassical approximation, spin transport properties can be described 
by the Boltzmann distributions of electrons and magnons [12]. The boundary conditions 
are to link the non-equilibrium electron distribution function of the metal to the magnon 
distribution function of the MIL. Within the model of the s-d exchange interaction (see 
Sec. Ill), the total angular momentum is conserved and thus for an ideal interface the first 
boundary condition would be 

Uo-)= 3t (o + ) (l) 

where j t is the total angular momentum current, and we assign the interface at x — 0. If we 
consider the left layer as a non-magnetic metal (x < 0), the angular momentum is carried by 
conduction electrons only and thus jt(0~) = ^(O - ) where j s denotes the conventional spin 
current density. For the MIL on the right, the angular momentum is carried by magnons only 
and thus it(0 + ) = J , m (0 + ) where j m corresponds to the magnon current density. Therefore, 
we may rewrite Eq. (1) as 

W-)= 3m (0 + ) (2) 

Note that for a magnetic metal, both spin and magnon current contribute to the total 
angular momentum current. 
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The second boundary condition is the relation among the electron spin accumulation 
5m 5 (0~), the magnon accumulation 5m m (0 + ), and the total angular momentum current 

it(o), 

G em Sm s {0-) - G me Sm m {0 + ) = j t (0) (3) 

where the two coefficients G em and G me will be calculated within the s-d model in the 
next section. The physics of this boundary condition is rather transparent: the first term 
represents the generation of the magnon current in the presence of electron spin accumulation 
and the second term describes the spin current produced by magnon accumulation. The 
combination of these two processes at the interface yields the total interface spin current. 
We immediately note that Eq. (3) is analogous to the case of spin current between two 
metallic layers in which the boundary condition is G cr /i cr (0 + ) — G^/i^O - ) = jo-(O) where \i G 
denotes the spin dependent chemical potential (a = ±1 or f (I) corresponds to spin-up 
(down)) which is proportional to spin accumulation, and G a characterizes the interfacial 
spin conductance [T3] . With this analogy, we may identify the coefficients G em (G me ) as the 
interface conductance for the conversion of the spin (magnon) accumulation to the magnon 
(spin) current; we simply call G em and G me spin convertance for convenience hereafter. 

We point out that the boundary condition, Eq. (3), is different from what we proposed 
in the earlier paper [12] where we related the spin and magnon accumulation via a local 
magnetic susceptibility. Clearly, such approximation corresponds to an ideal case in which 
the interface spin resistance is zero or the spin convertance is infinite. In the next section, 
we shall calculate these spin convertances and show that they are in fact finite and thus the 
magnon mediated electric drag effect predicted in Ref. [12] was overestimated by one order 
of magnitude. 



III. MICROSCOPIC CALCULATION OF SPIN CONVERTANCE 

We start with the s-d exchange coupling at a metal|MIL interface, 

H sd = -J sd J — Yl K c kt Ck/ ^ + a q cl; c kt)£k/= q +k (4) 

S k,q,k/ 

where cj^ (c^) and cj^ ( c ki) are the creation (annihilation) operators for spin-up and spin 
spin-down electrons respectively, (a q ) is the creation (annihilation) operator for magnons, 
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FIG. 1: Spin angular momentum transfer at a metal | MIL interface. The upper (lower) panel 
describes magnon current j e ^ m (spin current j m ^ e ) generated by spin accumulation 5m s (magnon 
accumulation 5m m ) at the interface. 

S is the spin per atom of the MIL, and N s is the number of atomic spins of the MIL at the 
interface. The exchange coupling strength J s d is given by the exchange integral with the 
overlapped wavefunctions of the conduction electrons and the magnetic ions. Since we do 
not know the detailed orbitals for the interface states, the magnitude of J s d at interface is 
less known compared to that in bulk materials and we will treat it as a parameter. 

The above exchange interaction gives rise to angular momentum transfer between the 
electron spins at the metallic side and the magnons at the MIL side. In equilibrium, the 
net spin current across the interface is zero. At non-equilibrium when there is a spin ac- 
cumulation at x = 0~ or a magnon accumulation at x = + , a net magnon/spin current 
may be present across the interface. In Fig. 1, we illustrate two angular momentum transfer 
processes. The total angular momentum current across the interface should be caused by 
both processes. We shall calculate them separately below. 

Magnon current generated by spin accumulation at the interface is defined as 



where <> refers to the thermal averaging over all states and Aj is the area of the interface. 
By explicitly working out the above commutator and by using the Fermi-golden rule, we 




have 



Je 



= 2 ^ (jf) 2 E [(^ + !)(! " Mfa ~ N ^ ~ Mht] <%k + E q - e*) (6) 

where 7V q and f^ a are the magnon and electron distribution functions respectively. We 
have considered a rough interface such that there is no correlation between the electron and 
magnon momenta for the magnon emission/absorption processes (i.e., we do not impose 
k' — k = q). We first consider the process in the upper panel of Fig. 1, i.e., magnon 
current due to electron spin accumulation. Accordingly, we take the equilibrium distribution 
function for magnons, i.e., 7V q = N® = [exp(£ , q //c j eT) — where the spin-wave energy 
is = Aq 2 + A p , the exchange stiffness is associated with the Curie temperature via 
A = 3kBT c al I /7r 2 (S + 1) [14J, and A^ is the spin wave gap due to magnetic anisotropy. 
The electron distribution function can be conveniently separated into equilibrium and non- 
equilibrium parts, 

h a = f° + ^[-S^(x)+g a (k,x)] (7) 

where is the Fermi distribution function, Sfi a (x) is the local variation of the chemical 
potential and <?o-(k, x) is the anisotropic part of the non-equilibrium distribution function 
(J d 3 kg a (k^x) = 0). By placing the above equilibrium magnon distribution function and 
non- equilibrium electron distribution function into Eq. (6), we arrive at 

je^m = G em Sm s (0~) (8) 

where we have defined Sm s = /x B ^ e (e j p)(5/x^ — 5fi\) as the spin accumulation with g e (sF) 
being the interface electron density of states at Fermi level. The spin convertance can be 
formulated by 



En 



it S f f E \ 

Gem = 2hk T Jsd9e(£F)a 2 Q M alj / dE^E^CSctf l^) (9) 

where aoM and a>oi are the lattice constants of the metal layer and the MIL respectively, 
g m (£q) is the interface magnon density of states, and E max (~ 3kBT c /(S + 1)) is the max- 
imum magnon energy If a parabolic magnon dispersion is assumed, then the dominant 
temperature dependence of G em is (T/T c ) 3 / 2 . The above result has already been obtained 

in dans]. 
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The spin current induced by magnon accumulation at the metal|MIL interface can be 
similarly calculated. We define this interface spin current as 

_ / 2 /% 



3 



m— Ye — 



ifiAi 



^a*a q , H sd 



(10) 



q 

After working out the ensured commutator, we find the spin current has exactly the same 
expression as Eq. (6); this is not surprising because the s-d interaction conserves the total 
angular momenta. To evaluate the spin current induced by magnon accumulation (see the 
process displayed in the lower panel of Fig. 1), we replace the electron distribution by the 
equilibrium value, f\^ a = /£, and separate the magnon density into equilibrium and non- 
equilibrium ingredients AT q = N® + SN^. We then find from Eq. (10), 

G me 5m m (0 + ) (11) 

where Sm m = (2/ig) J dE^g m {E^)8N^ is defined as the magnon accumulation. The spin 
convertance can be expressed as 

Gme = -^Jsd9l(^F)al M al I E m (12) 

with 

E-max 



J dEqg m (Eq)EqN® 

E m = ^ (13) 

J dE^E^N* 

A, 

where we have replaced the non-equilibrium magnon energy by the average magnon energy 
E m by assuming a near equilibrium magnon distribution. A rough estimation for simple 
parabolic bands of both magnons and electrons gives G me ~ (7rSal I /haoM)Jsd9e( e F) (^) 
where T F is the Fermi temperature of the metal layer. 

By combining Eq. (8) and (11), we attain Eq. (3) with the spin convert ances G em and 
G me given by Eqs. (9) and (12). 



IV. ROLE OF SPIN CONVERTANCE IN ELECTRICAL DRAG 

To experimentally realize the conversion between spin current and magnon current and 
to quantify the spin convertance, one needs to create a non-equilibrium condition such that 
a spin current or a magnon current can be generated, manipulated, and more critically, 
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FIG. 2: Schematics of three hypothetical devices: (a) NM|MIL|NM trilayers, (b) MM | MIL | MM 
trilayers, and (c) NM|MIL bilayers. In (a) and (b), a spin current is generated by an injected 
electric current via spin Hall effect, while in (c), a magnon current is induced by applying a 
thermal gradient. In all three cases, the magnetization directions of the MIL and the magnetic 
metal (MM) layers are oriented in +z. 

detected. The non-equilibrium states may be created in several ways. In this section, we 
study both electrical injection into a metal layer and thermal gradient across a MIL. In 
Fig. 2, we show three hypothetical devices to explicitly demonstrate the magnon-mediated 
electrical drag. 



A. NM|MIL|NM trilayers 

In Fig. 2(a), a magnetic insulator layer (MIL) is sandwiched between two non-magnetic 
metal (NM) layers. By applying an in-plane electrical current in the NM1 layer, a spin 
current flowing perpendicular to the layers would be generated due to the spin Hall effect. 
In this geometry, a partial spin current would flow into the MIL via transfer of spin current 
to magnon current. If the magnon diffusion length is larger than the thickness of the MIL, 
the magnon current would reach the other side of the MIL and subsequently, converts back 
to spin current in the NM2 layer. Finally, an electric current parallel to the layer is generated 
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owing to the inverse spin Hall effect [IT]. Such an electric drag phenomenon, namely, an 
electric current in one NM layer induces an electric current in the other when the two NM 
layers are separated by a MIL, would be a proof of the magnon/spin current conversion at 
magnetic interfaces. Although we have already calculated the drag coefficient in Ref. [12], 
we find the improved boundary conditions presented in this paper quantitatively modify the 
earlier result. 

Referring to the coordinate system in Fig. 2, one can establish the relation of the spin 
accumulation, spin current, magnon accumulation, and magnon current in each layer. For 
the NM1 layer with an applied in-plane current density we have 

5m s (x) = Ai exp(x/A s/ ) (14) 

where X s f is the spin diffusion length, A\ is a constant to be determined via boundary 
conditions. We have taken the thickness of the layer much larger than the spin diffusion 
length such that the term proportional to exp(— xj\ s f) has been dropped. The spin current 
flowing perpendicular to the plane of the layers is given by 

^) = -l sh ^-D s d -^ (15) 

where the first term represents the spin Hall effect: an electric current j^p in the ^-direction 
generates a transverse spin current proportional to the spin Hall angle ^y sh which is defined 
as the ratio of the spin Hall conductivity to the electric conductivity. Note that we have 
adopted e = \i B = 1 for notation convenience so that the electrical current and the spin 
current would have the same unit. The second term corresponds to the spin diffusion where 
D s is the spin diffusion coefficient which may be related to the conductivities = by the 
Einstein relation: = e 2 g e (sF)D s . For the MIL layer, we have 

Sm m (x) = A 2 exp(x// m ) + A 3 exp(-x// m ) (16) 

and 

dSiTi 

3m{x) = -D m m (17) 

where A 2 and A 3 are integral constants from the magnon diffusion equation, l m is the magnon 
diffusion length, D m is the magnon diffusion constant associated with the magnon diffusion 



length by l m = \J D m r t h with r t h being the magnon-nonconserving relaxation time [12] . For 
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the NM2 layer, we have 



5m s (x) = A4 exp(x/X s f) + A$ exp(— x/X 



sf) 



and 



j?\x) = -D. 



d8m s 
dx 



(18) 



(19) 



where A4 and A$ are two integration constants. The four boundary conditions of Eqs. (2) 
and (3) at the two interfaces x = and x — rf, along with the outer-boundary condition at 
x = d + L 2 where ji 2 \x = d + L 2 ) = 0, determine the five constants Ai (i = 1 — 5). After 
a straighforward algebra, we find the spin current density ji 2 \x) which in turn converts to 



an in-plane charge current in the NM2 layer via the inverse spin Hall effect, i.e., j 
lshJs 2 \x). Explicitly, 



X 



— a&sinh 


d-\-L2— x 

Xsf 


csch (ft) ^ 




sinh {£) + 


63 + a&coth (j^j 


cosh (£) 



(20) 



where we have introduced the dimensionless constants a = X s fG em / D s , b = l m G me /D m 
bi = 1 + a + 6 2 , 6 2 = a + a 2 , and 63 = (2 + a) 6. We may define an average electric current 
density by averaging over the thickness of the NM2 layer, = (I/L2) J dxji 2 \x). Then 



the ratio of the averaged current density to the injected current density, i.e., rj = 
can be obtained as 



77 



a&tanh(^) 



(21) 



L2 [&1 + k coth (^)] sinh (£) + [b 3 + ab coth (^)] cosh (^) 
The electrical drag coefficient 77 may be readily estimated. In the case of d <C Z m , 77 be- 
comes independent of G me , but increases with G em ; this is understandable since in this case 
the magnon current does not decay and thus the magnon accumulation is unimportant, rj 
depends predominantly on the efficiency of the magnon current generation by spin accu- 
mulation which is measured by G em . A quick numerical check also indicates that G em is 
usually larger than G me . We consider a trilayer of Ta|YIG|Ta whose material parameters 
at room temperature (T = 300K) are taken as follows: for the Ta layers [18], the conduc- 
tivity CTa = (190/ifl • cm) -1 , the spin diffusion length X s f = 5 nm and the spin Hall angle 
j sh = 0.15, the lattice constant a M = 3. 3^4, and the Fermi energy e F = 5 eV; for the YIG 
layer [19], the Curie temperature T c = 550 K ) the lattice constant a j = 12.376^4, the spin 
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FIG. 3: Electrical drag coefficient as a function of the NM2 (Ta) layer thickness for three different 
thicknesses of the MIL (YIG). See the main text for the parameters used in the figure. 

wave gap = 10 -6 eV : and the magnon relaxation time r t h — 10 -6 <s. In Fig. 3, we show 
7] as a function of the thicknesses of the NM2 layer for several different MIL thicknesses. 
Fig. 4 shows r] as a function of the interface exchange coupling J s d with several different 
magnon diffusion lengths. 

Finally, we discuss the sign of the drag current. The induced electric current always flows 
in the opposite direction of the injected electric current for any magnetization direction of 
the MIL . To see this, we first recall the spin Hall and inverse spin Hall effect in a single layer: 
an electric current induces a perpendicular spin current (spin Hall) which in turn produces 
an electric current (inverse spin Hall). The physical principle is that the combined actions 
of the spin Hall and the inverse spin Hall are to reduce the original driving current. Now 
consider the trilayer system. Since the spin current injected into the NM2 layer remains 
parallel to the spin current in the NM1 layer, the electric drag current in the NM2 layer 
must be antiparallel to the applied electric current in the NM1 layer. 

B. MM | MIL | MM trilayers 

Next, we consider a trilayer structure where the two metallic layers are magnetic, as shown 
in Fig. 2(b). Since the direct contact between the magnetic metal (MM) layer and the MIL 
would make it difficult to rotate the magnetization of each layer independently, one may 
insert a thin non-magnetic layer at the interface to break direct magnetic coupling. When 
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FIG. 4: Electrical drag coefficient as a function of J S( i/sf for three different magnon diffusion 
lengths of the MIL (YIG). See the main text for the parameters used in the figure. 

an in-plane current is applied to the MM1 layer, an anomalous Hall current perpendicular to 
the layers is generated if the magnetization is oriented in the z-axis. Although the physics 
of anomalous Hall and spin Hall effects are the same, the anomalous Hall current has both 
spin and charge currents. The charge current, however, is unable to penetrate the MIL; this 
leads to a charge accumulation at the interface so that the net charge current is exactly zero 
in the steady state. The spin current, on the other hand, is able to propagate into the MIL 
via the conversion to the magnon current, as discussed in the previous section. To gain a 
quantitative understanding, we carry out the following calculation. 

The x-components of spin and charge currents of the MM1 layer can be expressed as 

3s P s dx s dx lah3e [ } 

and 

Mi) _ D dSn™ dSmP (1) 

J e ~ » s Qx pV s Qx Vlah3e (26) 

where p = (c^ — c±) / (c^ + c±) is the spin polarization of the conductivity, ^ ah is the anomalous 
Hall angle defined as the ratio of the Hall conductivity to the longitudinal conductivity, 5n 
is the charge accumulation and is the current density applied in the ^/-direction as before. 
We have assumed a spin-independent spin diffusion coefficient D s . Since jl^ = 0, we may 
eliminate the charge accumulation term from Eq. (22) and get 

£<D = _(i _ p*)D s ^ - (1 - P W 1 ) (24) 
12 



For the MIL, Eqs. (16) and (17) remain valid, while for the MM2 layer, we similarly have 

fs m = -^-P 2 )Ds^- (25) 

By comparing Eqs. (24) and (25) with Eqs. (15) and (19), one should realize that the induced 
electric current in the MM2 layer can be simply obtained by replacing D s by (1 —p 2 )D s 
and j sh by (1 — p 2 )j a h m Eq. (20). Consequently, the electrical drag current is reduced by 
a fact of (1 — p 2 ) 2 for the MM | MIL | MM trilayer if one approximates ^y sh « 7 a/l . This might 
be counter-intuitive at first glance, since one would expect the magnetic metals to provide 
more spin signals. However, if we realize the interplay between the charge and spin currents, 
one can readily explain the above conclusion: consider the extreme case of p = 1, i.e., spin 
current generated by the anomalous Hall is fully polarized such that the spin current is 
same as the charge current. Since the charge current is completely blocked by the MIL, it 
is inevitable that the spin current is also being completely blocked. 



C. NM|MIL bilayers 

In this section, we consider a NM|MIL bilayer. In this case, the magnon current in 
the MIL is induced by a temperature gradient, see Fig. 2(c). From the magnon Boltzmann 
equation within the relaxation time approximation, the non-equilibrium magnon distribution 
is, 

XAT ^ dN * dT ^ 96 N « fort 

where denotes the x-component of the magnon velocity and r m is the magnon-conserving 
relaxation time. By defining the magnon current as j m = (2/i B ) J dqy^SNq and following 
the derivation in the Appendix A of the Supplemental Material of Ref . [12] , we find 

dT d5m m 
3m = ~^~ Dm ^h r (27) 

with 

3 
2 



V " " " (ir) (28) 



2^3(3 + lKr m fc|r c £ (T 
9irh 2 a 0I \T C 

and f = J °° dxx 3 / 2 e x /(e x - if ~ 3.4. The magnon accumulation satisfies the magnon 
diffusion equation whose solution can be taken as a simple form, 

Sm m (x) Bi exp(x// m ) (29) 
13 



where we have assumed the thickness of the MIL to be much larger than the magnon diffusion 
length l m and hence dropped the term exp(— x/l m ) in the solution. The spin accumulation 
in the NM layer can be written as 

5m s (x) = B 2 exp(x/X sf ) + B 3 exp(-x/A s/ ) (30) 

and the spin current perpendicular to the plane is given by j s = —D s d5m s /dx. These 
three integral constants Bi {i = 1,2,3) can be determined by the two interface boundary 
conditions Eqs. (2) and (3), along with the outer boundary condition j s (x = d + L 2 ) = 0. 
After a straightforward algebra, we get 

Js(x) = 7 \ ) 7 g (31) 

(1 + 6)sinh ( f^j + a cosh (^M ax 

Again, the above perpendicular-to-plane spin current can generate an in-plane electric 
current whose average density over the thickness of the NM layer can be obtained by 
j e = ( / y sh /L 2 ) J j s (x)dx. By taking the temperature gradient as a constant, we have 



Je 





cosh 


(ft) 


-1 


dT 


L 2 (1 + 6) s inh 


(ft) 


+ a cosh 


(ft) 


dx 



(32) 



The direction of the electric current is in the plane of the layer and perpendicular to the 
directions of the magnetization as well as the temperature gradient of the MIL. It is in- 
teresting to compare the current driven electric drag, Eq. (21), with the thermally driven 
electrical drag, Eq. (32). Firstly, in the former case, the electric drag is proportional to the 
square of the Hall angle because the first metallic layer converts the electric current to the 
spin current via spin Hall effect and the second metal layer converts the spin current into 
the electric current via inverse spin Hall effect, while in the latter case, the spin current 
is directly injected from the thermally driven magnon current and thus the drag current is 
linearly proportional to the spin Hall angle. Secondly, in the NM|MIL|NM case, both spin 
convert ances G em and G me are important, while for NM|MIL, the convertance relating the 
magnon accumulation to the spin current, G me , plays a dominant role. A rough estimation 
yields the induced current density in a Pt|YIG bilayer is about 10 A/ cm 2 for a moderately 
small temperature gradient of dT/dx = 10 K/cm if one chooses the following parameters: 
7 s/l 0.05 [20J, J s d 1 meV, L 2 X s f 7 nm, c Pt 0.1 (/jlQ • cm)' 1 , a 0M 3.9 A, 
a i = 12.376 i, r m = 10- 8 s [19J, r th = W~ 6 s,l m = l mm, S = § and T c = 550 K. 
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V. DISCUSSIONS AND SUMMARY 



We have investigated the spin transport across the interface between a metal layer and a 
MIL. The salient feature of our approach is that we have treated spin and magnon transport 
properties on an equal footing. Namely, the spin and magnon accumulations as well as 
the spin and magnon currents are described by semiclassical non-equilibrium distribution 
functions. In other approaches, for example, Xiao et al [211 E2] described the magnon 
density through a quasi-equilibrium effective magnon temperature which differs from the 
lattice temperature. Adachi et al [23l [24] considered the linear response theories and 
their numerical solutions [25] on the spin Seeback effect [26], [27] in ferromagnetic insulators. 
These approaches also provide alternative physical insights on the roles of magnons in non- 
equilibrium transport [28j . 

We thank S. Bender for pointing out an inconsistency in the approximations used in the 
derivation of the spin convertance G em in a previous version of the paper. This work is 
supported by NSF-ECCS. 
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